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Abstract. Families of exact solutions are found to a nonlinear modification of the Black- 
Scholes equation. This risk-adjusted pricing methodology model (RAPM) incorporates both 
transaction costs and the risk from a volatile portfolio. Using the Lie group analysis we 
obtain the Lie algebra admitted by the RAPM equation. It gives us the possibility to describe 
an optimal system of subalgebras and correspondingly the set of invariant solutions to the 
model. In this way we can describe the complete set of possible reductions of the nonlinear 
RAPM model. Reductions are given in the form of different second order ordinary differential 
equations. In all cases we provide solutions to these equations in an exact or parametric form. 
We discuss the properties of these reductions and the corresponding invariant solutions. 
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One of the most important problems at present is how to incorporate both the 
transaction costs and the risk from a volatile (unprotected) portfolio into the gov- 
erning Black-Sholes equation. In the pioneering work of Leland [ll|, devoted to the 
problem of option pricing in the presence of transaction costs, the idea of a periodic 
revision of a hedging portfolio was introduced. Leland assumed that the level of 
transaction costs is a constant, i.e. we have a market with proportional transaction 
costs. He reduced this problem to a nonlinear partial differential equation with an 
adjusted volatility. Leland claimed that the terminal value of the portfolio approx- 
imates the payoff as the length of a revision interval tends to zero. Later, Kabanov 
and Safarian [8] proved that Leland's conjecture based on approximate replication 
fails and his model has a non-trivial limiting hedging error relative to simulated 
marked prices (see as well the detailed discussion in [9|). Mathematical problems 
arise in the limiting cases as revisions become unboundedly frequent. As a practical 
matter, extremely frequent revisions will not be desirable and the average errors 
are less than one-half of one per cent of the price suggested by Leland's formula 
121 ] . Within this model Kratka [7j has suggested a mathematical method for pric- 
ing derivative securities in the presence of proportional transaction costs and he 
additionally took into account the risk of the unprotected portfolio in between the 
revisions. Jandacka and Sevcovic 0] modified Kratka's approach in order to derive 
a scale-invariant model. 
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In the model the risk from the volatile portfolio is described by the average value 
of the variance of the synthesized portfolio. The mathematical model was referred 
to as the risk-adjusted pricing methodology (RAPM) model. The RAPM model 
generalizes the famous Black-Scholes model for pricing of derivative securities. In 
the model setting both the transaction costs and the unprotected portfolio risk 
depend on the time interval between two transactions and minimizing of the total 
risk leads to the RAPM model. The model was studied recently with numerical 
methods in the case of European and American options fl~6l ] . 
We describe briefly the model settings. 

The authors assume that the stock price dynamics is given by the geometric 
Brownian motion 

S t = S e W ((p-a 2 /2)t + aW t ), (1) 

where {Wt,t > 0} is the Wiener process, p E R is the drift and a > is the 
instantaneous volatility of the asset, p, a are constants. It is assumed that the risk- 
free bond earns at a continuously compounded constant rate r. 
The time-steps At at which the portfolio can be hedged against the price change of 
the underlying asset St are non-infinitesimal and fixed. Additionally, the authors 
introduce the idea of a switching time t* for the last revision of the portfolio. This 
means that the time interval (0, T) is divided in two parts, in the first part (0,t*) 
the revisions of portfolio will be done regularly, and in the second one (t*,T) there 
are no revisions and correspondingly no transaction costs. It is assumed that the 
interval (t* , T) is very small and in this interval the price of the contingent claim 
u(S,t), t £ [t*,T] is defined as in the classical Black-Scholes formula (here T is 
the maturity time). It is assumed that the model (similar to Leland's model) does 
not include the cost of establishing the initial investor's portfolio composition. 
At time t the value of the dynamically hedged portfolio Vt is Vf = StSt + PtBt, 
where 5t is a number of units of the stock (a constant on each time interval At), 
B t is the value of the bond and fit is a number of units of the bond. We can put 
Bq = 1 without loss of generality and rewrite the previous relation in the form 
V t = 5tSt + /3te rt . The pair 4> = (5t,f3t) defines the self-financing hedging strategy 
that maintains the portfolio. 

The change of Vf in any time-step At is equal to AVf = V^_ At — V t = 
Pte rt (e rAt — 1) + 5t(St+At — St) — r^StAt. The total risk premium tr contains 
two parts rji = r^c + r VP- The transaction costs (TC) in this case are modeled by 
the expression 

CaS\uss\ „ S as k — Sud , nS 

rTC = ^r> c= ^^' (2) 



where C is the round trip transaction costs per unit dollar of transaction 11], 0], 
lQf ] and u(S, t) is the value function of the contingent claim with respect to the 
asset price S and time t. During the time-step At the portfolio is unprotected and 
the risk connected with a volatile portfolio (VP) is modeled by 

r VP = ^Ra 4 S 2 (u SS ) 2 At, (3) 

where R is a risk premium coefficient introduced in [7| and [6| and represents the 
marginal value of investor's exposure to a risk. The total risk premium depends on 
the time-lag At and it is a strong convex function between two consecutive port- 
folio revisions [la ]. To obtain a risk-adjusted Black-Scholes equation the authors 
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minimize the total risk premium m = r?c + r VP ■ They then obtain for the optimal 
time-lag the following value 



C 2/3 



a 2 {RV2^\Su ss \) 2 / 3 ' 



At pt 



Using Ito's formula the authors of [6] finally obtain the risk-adjusted pricing 
methodology model 



1/3 



u t + ^a 2 S 2 u S s(l ~ K Su ss)*) ~ ru + rSu s = 0, /i = 3 ( J , (4) 

where t G (0,t*) and the value t* is determined by the implicit equation T — t* = 
min5 > o At opt (S, t*). The equation represents a well-posed parabolic problem under 
the condition that 

Su ss (S,t)< (jX. (5) 

The condition © will not be fulfilled for usual Call and Put options at S = E 
and t —> T~ , where E is the strike price of the corresponding option. To avoid the 
singularities in the model the authors introduced the switching time t* such that 
condition ([5]) is satisfied by t = t*. The equation for t* which can be reduced to 
the form T — t* = CR~ 1 a~ 2 (for European Call and Put options) has a positive 
solution and the condition ([5]) is satisfied if 

^<cr 2 r, CR<*-. (6) 

From the analytical point of view this model is represented by a fully nonlinear 
parabolic differential equation (PDE) with a singular perturbation. Our goal is the 
study of the RAPM model with the methods of Geometrical Analysis. 



2. Symmetry properties 

Equation @ is the main subject of our investigations. The equation possesses a 
complicated analytical and algebraic structure. We provide the Lie group analysis 
of this equation with the goal of describing the complete set of symmetries of equa- 
tion (jl]) and to obtain possible reductions. Using the invariants of the subgroups 
of the symmetry group of the studied equation we reduce the partial differential 
equation to ordinary differential equations (ODEs). Solutions to these ODEs give 
us the invariant solutions to the nonlinear RAPM model in an analytical form. 
We obtain the symmetry group of the RAPM model in the way suggested by Sophus 
Lie and developed further in 14], [3] and 0]. We first find, using the Lie deter- 



mining equations, the Lie algebra L r of a dimension r admitted by the equation. 
Then we use an exponential map exp : L r — > G r and obtain the transformations of 
the symmetry grou p G r . To each subalgebra hi C L r corresponds a subgroup Hi 
of G r 0, H, Hi]. In most cases we do not need the explicit form of the group 



transformations and use directly the subalgebras hi of L r in order to reduce the 
RAPM model. 
In this way we prove the following theorem. 
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4 L. A. Bordag 

Theorem 2.1 : The equation ^ admits a four dimensional Lie algebra L4 with 
the following infinitesimal generators 

The commutator relations are 

[U U U 2 ] = -U 2 , [U 2 ,U 3 ] = ~rU 2 , (8) 

[U x ,U 3 ] = [U 1 , U4] = [U 2 ,U 4 ] = [U 3 ,U4} = 0- (9) 

The commutator relations ([8]) depend on the parameter r, i.e. on the interest rate 
included in the model. Depending on whether r = or r / 0, we obtain different 
commutation relations for the algebra generators of the Lie algebra L4. After the 
proper choice of generators we obtain, in both cases, isomorphic algebras. 
All four-dimensional real Lie algebras were classified by Patera and Winternitzs 



151 ]. We will use this classification and the corresponding notations for generators 
of L4. The algebra is spanned by the following generators L4 =< ei,e2,e3,e 4 >, 
which will have different meaning depending on the value of r. We denote a two 
dimensional Lie algebra spanned by two operators ei, e 2 with the unique non-trivial 
commutator [e\, e 2 ] = e 2 as L 2 . The algebra L 4 is a decomposable Lie algebra and 
can be written as a semi-direct sum 

L * = L 2(J)e 3 (J)e4, L 2 =< ei,e 2 >, [e 1 ,e 2 ] = e 2 . (10) 

Case r ^ 0. In the case r / the generators take the form 

.,_(,- m + U 3 = (r- 1)S± + (r - !)»!; + |, e 2 = U 2 = 

c a =rU 1 + U,=r S ± + ru± + ^, ^ = U t = S^. (11) 

Case r = 0. Using the previous notations we can represent L4 in the case r = 
in the form 

ei = -Ux = -S—-u—, e 2 = U 2 = — , 
ob ou ou 

«. = «=!, « = K. = S|;. (12) 



Patera and Winternitzs 15jJ looked for classifications of the sub-algebras into 
equivalence classes under their group of inner automorphisms. They also used the 
idea of normalization which guarantees that the constructed optimal system of 
subalgebras is unique up to the isomorphisms. 

This classification allows us to divide the invariant solutions into non-intersecting 
equivalence classes. In this way it is possible to find the complete set of essential 
different invariant solutions to the equation under consideration. We use this clas- 
sification and give a list of all non-conjugate one-, two- and three-dimensional 
subalgebras. The optimal normalized system of subalgebras to the algebra L4 is 
listed in Table [H 

In Tabled] we use the operators ei,e 2 ,es given by (fl2|) if r = and by ([TT]) if 
r / 0. 
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Tabic 1. I X 511 The optimal system of subalgcbras hi of the algebra L4 where a £ E, e — 

±1, <t> g [0,71-]. 



Dimension 


Subalgebras 


1 


h-\ =^ po /) =^ po cos { th} -\- p /i { fh) ~> 
/13 =< ei + a(e3 cos (</>) + e4 sin (0)) >, 
hn=< e2 + e(e3 cos (<i) + e4 sin (<i>)) > 


2 


h 5 =< ei + a(e3 cos (</>) + e4 sin (</>)), e2 >, /i6=< e 3, e 4 >, 
hi =< ei + a(e3 cos (<^)) + e4 sin (</>)), e3 sin (0) — e4 cos >, 
^8 =< e 2 + f (^3 cos (</>) + e4 sin (</>)), e3 sin (0) — cos (</>) >, 
hg =< e2, e3 sin (tj>) — pa cos (0) > 


3 


h\o =< ei,e3,e4>, hn =< e?.,Pi,PA >, 
hi2 =< ei + a(e3 cos (0) + e4 sin (</>)), e3 sin (0) — e4 cos (</)), e2 > 



In correspondence with the set of subalgebras listed in Table [H we obtain the 
complete set of invariant functions and reduce equation (J5|) to different ODEs using 
these functions as dependent and independent variables. 



3. Group-invariant reductions provided by the one-dimensional symmetry 
subgroups in the case r 7^ 

In this section we study the symmetry reductions of the RAPM model which 
we obtain using one of the one-dimensional symmetry subgroups Hi,i = 1,...,4. 
These symmetry subgroups Hi C G4 are generated by the corresponding 
subalgebras hi,i = 1,...,4 listed in Table Q] by a usual exponential map. We 
skip the study of invariant reductions to the two and three dimensional sub- 
groups listed in Table Q] because they only give trivial results for the RAPM model. 



Case H\. This one-dimensional subgroup H\ is generated by the subalgebra 

7 rt ^ 

hi =< e 2 >=< e rt — > . 

ou 

It describes a gauge (or evolutionary) symmetry of the equation. It means that to 
each solution to equation (JJJ we can add a term ae rt , where a is arbitrary constant. 
The new function u(t, S) — > u(t, S) + ae rt is then still a solution to the equation. 
This symmetry does not give rise to any invariant reductions of equation (J3J). 

Case H 2 . We look for the invariants of the subalgebra h 2 =< e^cos(cf)) + 
e4 sin {(f)) >. In the variables (t, S, u) we obtain that h 2 has the form 

d d d 

h 2 =< cos(<ft)— +rcos((ft)S— + (cos((/)) r u + sin(<ft)S)— >. (13) 

The invariants z, w of the corresponding subgroup H 2 C G4 can be chosen in the 
form 

z = Se~ rt , w = ---lnS, r / 0, r = tan U), cj> € [0, tt], 6 4 tt/2. (14) 
S r 



We take the invariants z, w as the new independent and dependent variables, re- 
spectively, then the PDE (J5]) is reduced to the ordinary differential equation of the 
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following form 

(r + rz(zw zz + 2w z )) (l - fir~3 (r + rz(zw zz + 2w z ))^ + ^ = 0, (15) 

r / 0, r = tan ((f)), G [0,7r], (j> ^ 7r/2. 

This second order differential equation can be reduced to a first order equation by 
the substitution w z {z) = v(z) which has the form 

l-tir—t(T + r(z 2 v) z )~ 3 )+-^ r = 0. (16) 

From this equation it follows that the expression (z 2 v) z is a constant. If we denote 
(r + r(z 2 v) z ) 1 ^ 3 = p(z), then for the value p(z) we obtain an algebraic equation of 
the fourth order 

p 3 (l - fir-^p) + -JJ = 0. (17) 

This equation has four roots qi,i = 1, ... ,4. In dependence on the values of the 
constants /u and r some of these roots are real. We denote the real roots by ki. 
To find solutions to the ODE (fT5|) we have just to integrate two simple first order 
differential equations 

r + r(z 2 v) z = kf, w z {z)=v{z). (18) 

Then to each root ki the corresponding solutions to equation (|15p are given as two 
parametric families of functions 

k 3 

u(S, t) = ^SlnS- (kf - r)tS + aS + c 2 e rt , (19) 
r 

where ci, C2 GK,r/0,r = tan ((f)), (f) £ [0, 7r], <f) / ir/2. 

Case H%. The subalgebra /13 is spanned by the generator e\ + a(e^cos((j)) + 
e4 sin ((f))). In the variables (t, S, u) it means that we have to do with the subalgebra 
of the form 

(9 d 
h 3 =< (l + acos(0))— + ((r-l) + arcos + (20) 

d 

((r — l)u + a(cos ((f>)ru + sin ((f>)S)— > . 

ou 

The two first invariants of the corresponding subgroup H% are given by z, w which 
are connected to variables (t, S, u) by 

z = Se- (r+ ^ )t , u(S,t) = Sw(z) +(SlogS, (21) 

where the constants are 7 = (1 + acos(0)) _1 , C = r(i+acos(4>))-i > Q ^ & ^ Pi 71 "]- 
Using these expressions we reduce the RAPM equation to an ordinary differential 
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equation of the form 

2 

y (z(zw) zz + (l - n (z(zw) zz + C)")+rC~ jzw z = 0. (22) 

The solutions to this equation can be given in the parametric form 

d9 



z(9)=exp (J 



"W-fwf^-c (23) 

where 9 £ R is a parameter and qt{9) is one of the real roots of the fourth order 
algebraic equation 

2 

^-k i {ef{i-^k i {e)) + rc- 1 e = ^). (24) 



Case H4. The subalgebra hi is spanned by the generator e<i + a(e3cos(0) + 
e4 sin (</>)). In terms of the variables (t, S, u) it means that we are dealing with the 
subalgebra of the form 

d d d 

hi=< ecos((p)— + ercos (</>)S— + (e rt + e(cos (^)ru + sin (0)5)— >. (25) 
at ob au 

The invariants of the corresponding subgroup H4 are z and w, where 



z = Se~ rt , u(S,t) = Sw(z) + ( ~ + ^— z~A Slogff, (26) 



r r cos( 



with t = tan(^), 4> G [0, 7r], 7^ 7r/2 and e = ±1. We take these invariants as new 
invariant variables and reduce equation §S§ to an ODE of the following form 

c 2 //x t e \ A / , . t e 



— z(zw) zz H 1 — - 1 - A* z(ot) 22 H h 



2 \ r rzcos(0) / V \ r rzcos(0) 

+r + ^— = 0. (27) 

Z COS(0) 



If we denote p(z) = yz(zw) zz + 7 + rz C o S (^) ) then for the value we obtain 
an algebraic equation of the fourth order 

p 3 (z) (1 - w(z)) + %+ 2 26 = 0. (28) 
cr z cos(0J 

This equation has four roots which we denote % , i = 1, . . . , 4 as in the case i?2 • 

Remark. The roots in this equation differ from the roots of equation (I24j) 
or (I17p , Still, we denote here (and later) all real roots of a fourth order algebraic 
equation by ki to show the similar structure of solutions. 
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Then to each root ki(z) the corresponding solutions to equation ([4]) are given as 
two-parametric families of functions 

\3 



u(S.I) =e rt I ( / ^-d: j d: + S (rl + r, ) 



Id 

eH ( ^TJ^t + c 2 ) , (29) 



v COS I 

where r = tan(</>), z = Se~ rt , (j) E [0, ir], <p ^ ir/2, c\,C2 S R and e = ±1. 
The special case of invariant solutions. 

In some cases it is more rewarding not to take one of the classical representatives 
listed in Table Q] of the non-conjugated subalgebras but rather turn to an equivalent 
one which gives us a simpler ODE. Let us take a one-dimensional subalgebra of 
the form h =< e\ + ae<i >, where e%, ei are defined by (jlip . The invariants of the 
corresponding subgroup H are defined by the infinitesimal generator 

U = e 1 + ae 2 = (r-l)U 1 + U 3 + aU 2 , (30) 

and can be chosen in the form 

z = ^ e ~('- 1 )*, w = u(S,t)e~ (r ~ 1)t -ae l . (31) 

Remark. In the case r = 1 the dependence of the invariants z,w on t will be 
trivial. It means then that z = S is an invariant and w = u + ae . On the other 
hand, the value r = 1 implies that on the market 100 per cent interest rates are 
accepted. This is certainly a case which can not be modeled with the RAPM 
model. We can, therefore, exclude the case r = 1. 

We use these invariant functions z and w to reduce the original equation Q to 
the ODE of the form 

-w + zw z + ^a 2 z 2 w zz (1 - fi(zw zz ) 1/3 ) = 0. (32) 

It is easy to see that this equation does not depend on the arbitrary parameter a 
which is included in ([3D]) . The second order ODE (f52"j) can be reduced to a first 
order one 

vn-l«^ = ~ (33) 
a z z 

by the substitution 

v(z,w) = zw z — w. (34) 

Equation (|33p has a parametric solution. We obtain this solution in the following 
way. We rewrite equation (|33|) in the form 

2 

v(z) = - a —z (v z - fiv 4 J 3 ) = G(z,v z ), (35) 
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then the parametric solution to this equation is given by the solution to the system 
of equations 

v(e) = G(z(e),e), ^ = e G ffJ e) =~ , ^V/^" ' (36) 



where 6 £ M. is a parameter. The system (|36|) and correspondingly equation ()33[) 
have the following solution 

v (9) = ~z{O){e-y0V% z(e) = Cl (i-^ (i-f^elfj + "\~^\ (37) 

where 7 = (l — and c\ = const. Using the parametric solution (p7|) to (f33j) 



2 r 

we obtain the parametric solution to (|32|) . We used the substitution (j34"j) which 
now takes the form 

v (0) = z(6)w z - w = (In z(6>))0 1 w e - w. (38) 

This is a linear first order differential equation for the function w(t) and together 
with the parametric representation of z{6) (|37p the solution to this equation gives 
us the parametric solution to (|32l) 



w{6) = z{6) (ca + g(0)) , c 2 = const, (39) 
where the function g{9) is given by 

2 fi \ / 2\ — 3 / 2\ / 2\2 / 2 



Expressions (|39~]) and ([37]) give a parametric representation of a solution w(z) to 
equation (f32j) . 



4. Group-invariant reductions provided by one-dimensional symmetry 
subgroups in the case r — 

We repeat the procedure of constructing the invariant solutions to the RAPM 
model in the case r = 0. The general structure of the optimal system of sub- 
algebras is the same in both cases but the form of infinitesimal generators differ. 
The invariants and the reductions therefore take another forms. 

Case H® . The generator of the subalgebra has a very simple form e 2 = ^7 in 
the case r = 0. This means that we are dealing with a subgroup of translations in 
the ^-direction. Hence, to each solution to equation (|4|) with r = 0, we can add an 
arbitrary constant without destroying the property of the function to be a solution. 
This subgroup does not provide any reduction. 
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Case H® • The subalgebra h® has the form h® =< e.3 cos (<f>) + sin (</>) >, it 
means that in terms of the variables (t, S, u) we have the subalgebra of the following 
type 

h° 2 =< cos (</»)- + 8^(0)5— >. (40) 

The invariants of the subgroup H® are given by 

Z = S, w = u(S,t) -rtS, t = tan (<£), <j) G [0, vr], / vr/2. (41) 

If we use the variables z, w as new independent and dependent variables we obtain 
the following reduction of the RAPM model (J4]) with r = 



2 



^(1-^(2^+7 = 0, r = tan (0), G [0, vr], <f> + vr/2. (42) 



We denote {z w zz ) 3 = p{z) and obtain for the value p{z) an algebraic fourth order 
equation 

p 3 (1 - up) + % = 0. (43) 

As before we denote the real roots of this equation by ki. To find solutions to the 
ODE (I42D we have just to integrate twice 

z w zz = kf. (44) 

Then the corresponding solutions to equation (|42p are given by 

u(S,t) = kfS(\nS-l)+ TtS + Cl S + c 2 , (45) 

where r = tan((j)), c\, cq G M, 4> € [0, 7r], 7^ 7r/2. 

Case i/g. The subalgebra /ig for r = has the form 

h° 3 =< acos(0)--5— + ( as in(0)S-u) — >, (46) 

where a G i2, </> G [0, 7r] are parameters. The invariants z, w of the group H® are 
given by the expressions 

z = Se 5t , u(S,t) = Sw(z) + (S log S, 

where the parameters are defined as 

S = (a cos^y 1 , ( = asm(4>), a G R, a / 0, (j) G [0,vr], <f> / vr/2, (47) 

and the reduced equation takes the form 
2 

y (z(zw) M + C) (l - M («H« + 0*) + 5zw z = 0. (48) 
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The solutions to this equation can be represented in the parametric form (123|) . 
where k{(v) is one of the real roots of the equation 

2 

^k i (vf(l-Lik i (v)) + 5v = 0, (49) 

and the parameter 5 is defined in (|47p . 

Case H®. The subalgebra h® for r = has the form 

/i°=< e cos (</>)— + (l + esin(0)S)^ >, (50) 

where e = ±1, G [0, tt] are parameters. 
The invariants z,w of this subgroup i/4 are given by the expressions 

z = S, w(z) = u(S,t) — rtS ~T7T ? t = tan (</>), (51) 

COS (<?)) 

and the RAPM model is reduced to the ODE of the form 

2 

?-z 2 w zz (l- n(zw zz )i) +tz+ — ^ = 0, (52) 
2 V / cos 



where r = tan(^), <f> G [0, 7r], 7^ vr/2, e = ±1. The structure of equation (i52j) 
is very similar to previous cases and we can use similar tools to solve it. We first 
substitute (zu; 22 ) 1//3 = p(z). Then for the function p(z) we obtain a fourth order 
algebraic equation but now its coefficients depend on the variable z 

p(zf (1 - np{z)) + % + 2 2e =0, (53) 
cr z za A cos (0) 

where r = tan (</!>), G [0, 7r], <fi ^ tt/2, e = ±1. For each real root ki(z) of this 
equation we have then to solve a linear ODE 

z w zz = h{zf. (54) 

The corresponding invariant solutions to (jU) then have the form 

u(t, S) = J ^-ds\ dS + tan (<f>)tS+ + ciS + c 2 , (55) 

where ci, c 2 £l, G [0,7r], ^ 7r/2, e = ±1. 

The expressions for these solutions are rather lengthy and because of which they 
are omitted here. 



5. Conclusion 

In the previous sections we found the complete series of reductions of the RAPM 
model. In this way the partial differential equation (j4j) is reduced to ordinary dif- 
ferential equations. Using the optimal system of subalgebras (Table [T]) allowed us 
to present the complete set of the non-equivalent reductions of equation (j3|) up 
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to the transformations of the group G4. In all cases it is possible to solve these 
ODEs and to obtain the exact or parametric representations of solutions to the 
RAPM model. We deal with the very seldom case that we can compare structures 
of non-equivalent invariant solutions since they are given in exact or parametric 
forms. Each of these solutions contains two integration parameters and some free 
parameters connected with the corresponding subgroup. This reasonable set of pa- 
rameters allowes one to approximate a wide class of boundary conditions. 
The RAPM model ([4]) possesses a non-trivial analytical and singular-perturbed al- 
gebraic structure. There exist rather few methods to study equations of such high 
complexity. An application of both analytical and numerical methods to singular- 
perturbed equations is a highly non-trivial task. The RAPM model was studied 
before in detail with numerical methods in 0] and in (l6l | . The authors of 0] derive a 
robust numerical scheme for solving equation and perform extensive numerical 
testing of the model and compare the results to real market data. In fl~6l ] Sevcovic 
studies the free boundary problem for the RAPM model and provides a descrip- 
tion of the early exercise boundary for American style Call options with floating 
strike. He proposed a numerical method based on the finite difference approxima- 
tion combined with an operator splitting technique for numerical approximation of 
the solution and computation of the free boundary condition position. 
On the other hand the Lie group analysis of the RAPM model which we provide 
in this paper gives us a more general, alternative point of view on the structure of 
this equation. It opens the possibility to exploit the Lie algebraic structure of the 
equation and may be helpful to improve another methods. 
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